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ABSTRACT 

We aim at studying the gamma-ray burst GRB 081008 environment by analysing 
the spectra of its optical afterglow. UVES/VLT high resolution spectroscopy of 
GRB 081008 was secured ~ 5 hr after the Swift-BAT trigger. Our dataset comprises 
also three VLT/FORS2 nearly simultaneous spectra of the same source. The avail- 
ability of nearly simultaneous high and low resolution spectra for a GRB afterglow 
is an extremely rare event. The GRB-Damped Lyman Alpha system at z = 1.9683 
shows that the interstellar medium (ISM) of the host galaxy is constituted by at least 
three components which contribute to the line profiles. Component I is the redmost 
one, and is 20 km/s and 78 km/s redward component II and III, respectively. We 
detect several ground state and excited absorption features in components I and II. 
These features have been used to compute the distances between the GRB and the 
absorbers. Component I is found to be 52 ± 6 pc away from the GRB, while com- 
ponent II presents few excited transitions and its distance is 200lgQ pc. Component 
III only features a few, low ionization and saturated lines suggesting that it is even 
farther from the GRB. Component I represents the closest absorber ever detected 
near a GRB. This (relatively) low distance can possibly be a consequence of a dense 
GRB environment, which prevents the GRB prompt /afterglow emission to strongly 
affect the ISM up to higher distances. The hydrogen column density associated to 
GRB 081008 is logA^/cnr 2 = 21.11 ± 0.10, and the metallicity of the host galaxy is 
in the range [X/H] = -1.29 to -0.52. In particular, we found [Fe/H]= -1.19±0.11 and 
[Zn/H]= —0.52 ±0.11 with respect to solar values. This discrepancy can be explained 
by the presence of dust in the GRB ISM, given the opposite refractory properties of 
iron and zinc. By deriving the depletion pattern for GRB 081008, we find the optical 
extinction in the visual band to be Ay ~ 0.19 mag. The Curve of Growth analysis 
applied to the FORS2 spectra brings column densities consistent at the 3a level to 
that evaluated from the UVES data using the line fitting procedure. This reflects the 
low saturation of the detected GRB 081008 absorption features. 

Key words: gamma-rays: bursts - ISM: abundances - line: profiles - atomic data. 



1 INTRODUCTION 

Long gamma-ray bursts (GRBs) are powerful explosions oc- 
curring at cosmological distances, linked to the death of 
massive stars. The gamma-ray (or prompt) event is followed 
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by an afterglow at longer wavelengths, which is crucial in 
order to understand the physics of these sources, but also 
to investigate the nature of the interstellar medium (ISM) 
of high redshift galaxies. Before GRBs, such studies made 
use of Lyman-break galaxies (LBGs, see e.g. Steidel et al. 
1999) and galaxies that happen to be along the lines of 
sight to bright background quasars, commonly referred to 
as QSO-Damped Lyman Alpha (DLA) systems. However, 
both classes are entangled by selection effects. In fact, LBGs 
fall in the bright end of the galaxy luminosity function and 
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Table 1. UVES and FORS2 setups. 



Instrument 


Setup (nm) 


Time from burst (hr) 


Exposure (s) 


Wavelenj 


5th (A) 


Slit width 


Resolution 


s/r 


s 


UVES 


Die 1, 346 


4.30 


1800 


3300 - 


3870 


1" 


40 000 


~ 3 - 


- 5 


UVES 


Die 1, 580 


4.30 


1800 


4780 - 


6810 


1" 


40 000 


~ 5 - 


- 8 


FORS2 


600B+22(A) 


4.37 


900 


3300 - 


6300 


1" 


780 


~ 35 - 


- 50 


FORS2 


600B+22(B) 


4.63 


900 


3300 - 


6300 


1" 


780 


~ 35 - 


- 50 


FORS2 


600B+22(C) 


4.88 


900 


3300 - 


6300 


1" 


780 


~ 35 - 


- 50 


FORS2 


A+B+C 


4.63 


2700 


3300 - 


6300 


1" 


780 


~ 60 - 


- 80 



may not entirely represent typical high-redshift galaxies. On 
the other hand, QSO sightlines preferentially probe galaxy 
halos, rather than bulges or discs, for cross-section effects 
(Fynbo et al. 2008). Indeed, Savaglio et al. (2004; 2005) stud- 
ied the ISM of a sample of faint AT-band selected galaxies 
at 1.4 < 2 < 2.0, finding Mgll and Fell abundances much 
higher than in QSO systems but similar to those in gamma- 
ray burst hosts. Unfortunately, these galaxies are too faint 
to be spectroscopically studied up to higher redshifts, us- 
ing 8m class telescopes. In this context, long GRBs can be 
used as torchlights to illuminate the high-redshift ISM, and 
thus represent an independent tool to study high-redshift 
galaxies. 

Several papers report a metallic content in GRB host 
galaxies in the range 10~ 2 — 1 with respect to solar val- 
ues (see e.g., Fynbo et al. 2006; Savaglio 2006; Prochaska 
et al. 2007). The GRB host metallicity is thus on average 
higher than in QSO-DLA systems, supporting the notion 
that GRBs explode well within their hosts. Since long GRBs 
are linked to the death of massive stars, they are though to 
originate in molecular clouds. In this scenario, absorption 
from ground-state and vibrationally excited levels of H2 and 
other molecules is expected, but not observed (Vreeswijk 
et al. 2004; Tumlinson et al. 2007). The non-detection of 
these molecular states (with the exception of GRB 080607, 
see Prochaska et al. 2009; Sheffer et al. 2009) could be a 
consequence of the intense UV flux from the GRB afterglow, 
which photo-dissociates the molecules. However, molecular 
hydrogen is not detected in QSO-DLA either (e.g., Noter- 
daeme et al. 2008; Tumlinson et al. 2007), possibly indicating 
that these molecules are just hard to see at high redshift. 

This is just an example of how a GRB can modify its 
surrounding medium. The most impressive manifestation of 
the transient nature of GRBs in optical spectroscopy is the 
detection of strong absorption features related to the ex- 
cited levels of the OI, Fell, Nill, Sill and CII species and 
their time variability (Vreeswijk et al. 2007). This variability 
can not be explained assuming infrared excitation or colli- 
sional processes (Prochaska, Chen & Bloom 2006; Vreeswijk 
et al. 2007; D'Elia et al. 2009a), thus excitation by the in- 
tense GRB UV flux is the leading mechanism to produce 
these features. In this framework, the GRB/absorber dis- 
tance can be evaluated comparing the observed ground state 
and excited level abundances with that predicted by time- 
dependent photo-excitation codes. This distance turns out 
to be in the range ~ 0.1 — 1 kpc (Vreeswijk et al. 2007; 
D'Elia et al. 2009a,b; Ledoux et al. 2009). 



Within the described framework, the best and most 
complete tool to perform these kind of studies is high reso- 
lution spectroscopy. In fact, it is the only way to disentangle 
the GRB interstellar medium in components and to separate 
the contribution to the absorption coming from the excited 
levels from the ground state ones. In addition, a high spec- 
tral resolution allows us to check for saturation of lines a 
few km s _1 wide that may appear unsaturated in lower res- 
olution spectra (see e.g. Penprase et al. 2010). 

In this paper we present data on GRB 081008, observed 
both in high and low resolution UVES and FORS2 at the 
VLT. The paper is organized as follows. Section 2 summa- 
rizes the GRB 081008 detection and observations from the 
literature; Sect. 3 presents the UVES observations and data 
reduction; Sect. 4 is devoted to the study of the features from 
the host galaxy, in particular their metallicity and distance 
from the GRB explosion site; Sect. 5 presents the FORS2 
data and makes a comparison with the UVES ones; finally in 
Sect. 6 the results are discussed and conclusions are drawn. 
We assume a concordance cosmology with Ho = 70 km s _1 
Mpc -1 , ttm = 0.3, Q A = 0.7. Hereafter, with [X/H] we refer 
to the X element abundance relative to solar values. 



2 GRB 081008 

GRB 081008 was discovered by Swift/BAT on October 8, 
2008, at 19:58:29 UT, and was detected by both the XRT 
and the UVOT instruments (Racusin et al. 2008). The 
UVOT magnitude in the white filter was reported to be 
15.0 at 96 s from the trigger. The afterglow was also de- 
tected in all filters (from B to K) by SMARTS / ANDIC AM 
~ 4 hr post burst (Cobb 2008). The redshift was secured 
by the Gemini-South/GMOS, which observed the after- 
glow 5 hr after the Swift trigger, reporting a redshift of 
z — 1.967 (Cucchiara et al. 2008a). This value was later 
confirmed by our VLT/UVES+FORS2 data (D'Avanzo et 
al. 2008). The host galaxy was identified in the Gemini- 
South/GMOS acquisition image, and spectroscopically con- 
firmed to be at the GRB redshift. The host of GRB 081008 
has R = 20.75 ± 0.01, which corresponds to an absolute 
AB magnitude of —21.5 (Cucchiara et al. 2008b). A multi- 
wavelength study of the prompt event and the early after- 
glow phase of GRB 081008 is reported by Yuan et al. (2010, 
hereafter Y10), which present Swift (BAT+XRT+UVOT), 
ROTSE-III and GROND data. 
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Table 2. Rest frame equivalent widths of the UVES features. 



Species 


Transition 


w r (A) 


A W r (A, Ict) 


OI 3 P 2 (g.s) 


1302 


0.57 


0.05 


Allien (es) 


1670 


0.67 


0.01 


AlIII 2 S 1/2 (g.s) 


1854 


0.22 


0.01 


1862 


0.13 


0.02 


SiII 2 P« /2 (g.s) 


1260 


0.63 


0.07 


1 808 


o 9n 
u.zu 


O 09 


bUi ^3/2 (i > 


1264 


0.63 


0.07 




1816 


0.06 


0.02 


CrIP5i /2 (g-s.) 


2056 


0.20 


0.02 


2062 


0.14 


0.02 




2066 


0.11 


0.02 


FeIIo°D 9/2 (g.s) 


2249 


0.12 


0.01 


2260 


0.20 


0.02 


FcIIa fi D 7/2 (1*) 


1618 


0.05 


0.01 


1621 


0.11 


0.01 


FcIIa 6 Dr/, (2*) 


1629 


0.03 


0.01 


FcIIa D D 3 / 2 (3*) 


1634 


0.03 


0.01 


1636 


0.03 


0.01 


lellosa -TQ/2 [o ) 


1637 


U.04 


0.01 


1612 


0.12 


0.01 




1702 


0.21 


0.02 


FeIIa 4 D 7/2 (9*) 


1635 


0.03 


0.01 


NiII 2 D 5/2 (g.s.) 


1741 


0.07 


0.02 


NiII 4 F 9/2 (2*) 


2166 


0.19 


0.01 


2217 


0.27 


0.01 




2223 


0.09 


0.02 


ZnII 2 5 1/2 (g.s.) 


2026 


0.19 


0.02 


2062 


0.09 


0.02 



3 UVES OBSERVATIONS AND DATA 
REDUCTION 

The GRB 081008 afterglow was observed with the high res- 
olution UV-visual echelle spectrograph (UVES, Dekker et 
al. 2000), mounted at the VLT-UT2 telescope, in the frame- 
work of the ESO program 082.A-0755. Observations began 
on the 9 th October 2008 at 00:16:43 UT (~ 4.25 hr after the 
Swtft/BAT trigger), when the magnitude of the afterglow 
was R ~ 18.5. Data were acquired under good observing con- 
ditions, with seeing ~ 0.7. Only the UVES-dichroic-1 (red 
and blue arm) was used due to observational and scheduling 
constraints. The net exposure time of the observation is 30 
minutes. The slit width was set to be 1" (corresponding to 
a resolution of R = 40000) and the read-out mode was re- 
binned to 2 x 2 pixels. The spectral range of our observation 
is ~3300A to ~3870A, ~4780A to ~5750A, and ~5830A to 
~6810A. Table 1 makes a summary of our observations. 

The data reduction was performed using the UVES 
pipeline (version 2.9.7, Ballester et al. 2000). The signal- 
to-noise ratio per pixel is ~ 3 — 5 in the blue arm and 



~ 5 — 8 in the red one. The noise spectrum, used to de- 
termine the errors in the best-fit line parameters, was calcu- 
lated from the real, background-subtracted spectrum, using 
line-free regions to evaluate the standard deviation of con- 
tinuum pixels. Since the noise spectrum has been produced 
after the pipeline processing and the background subtrac- 
tion, it takes into account possible systematic errors com- 
ing from the data reduction process. Fig. 1 shows the full, 
smoothed and normalized UVES spectrum. 



4 UVES DATA ANALYSIS 

The gas residing in the GRB host galaxy is responsible 
for many features observed in the GRB 081008 afterglow 
spectrum. Metallic features are apparent from neutral (OI) 
and low-ionization (A1II, A1III, Sill, CrII, Fell, Nill, Znll) 
species. In addition, strong absorption lines from the fine 
structure levels of Sill, Fell and from the metastable lev- 
els of Fell and Nill are identified, suggesting that the in- 
tense radiation field from the GRB excites such features. 
Table 2 gives a summary of all the absorption lines due to 
the host galaxy gas and report their rest frame equivalent 
widths (W r ). The spectral features were analyzed with FIT- 
LYMAN (Fontana & Ballester 1995), using the atomic pa- 
rameters given in Morton (2003). The probed ISM of the 
host galaxy is resolved into two main components separated 
by 20 km s -1 (Figs. 2 and 3). The wealth of metal- line tran- 
sitions allows us to precisely determine the redshift of the 
GRB host galaxy. This yields a vacuum-heliocentric value of 
z — 1.9683 ± 0.0001, setting the reference point to the red- 
most component (hereafter component I). The absorption 
features have been fitted with Voigt profiles, fixing the red- 
shift of the two components when studying different lines. 
All transitions appear to be nicely lined-up in redshift, with 
the exception of component II of SiIIA1808. We attribute 
this misalignment to a contamination of another feature, 
and fit just the redmost side of component II. All ground 
state and metastable species present absorption features in 
both components, while fine structure levels in component 

I only. Two sharp features can be seen at v = ±80 km s _1 
from the Fell a 6 D 5 / 2 line (Fig. 2). They can not be sepa- 
rated in the FORS2 spectrum, thus we can not safely as- 
sess if they are real or not. The Doppler b parameter has 
been linked between different excited transitions belonging 
to the same species. A small variation of the Doppler pa- 
rameter is allowed among different species, but the fits are 
quite good even fixing it. The values for components I and 

II are ~ 10 and ~ 20 km s _1 , respectively. An exception 
to this behaviour is represented by component II of Znll. 
In order to obtain a good fit, a b ~ 4 and ~ 50 km s is 
required for component I and II, respectively. This large b 
value in component II is necessary to adequately fit what 
appears as a low level of the continuum in particular in the 
ZnIIA2026 feature. The large difference between the b 
parameters deduced for Znll and CrII is odd, but 
we do not have a simple explanation for it. The col- 
umn densities and b parameters for all the elements and ions 
of the host galaxy's absorbing gas are reported in Table 3. 
A third component is actually identified at —78 km s _1 for 
some low-ionization lines only, i.e., the OIA1302, A1IIA1670 
and SiIIA1260, and for the fine structure level SiIIA1264 (see 
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Figure 1. The full, smoothed and normalized UVES spectrum. Solid lines indicate the noise level as a function of wavelength. 



6800 



Fig. 4). These lines are however heavily saturated, and their 
column densities reported in the table just set a lower limit 
to the true values. The reported upper limits are at the 90% 
confidence level. 



4.1 Abundances 

The GRB 081008 redshift was high enough to allow the 
hydrogen Lya line to enter the UVES spectral window. 
Unfortunately, the UVES spectrum is extremely noisy in 
this region, and the derived hydrogen column density, 
log(7VHi/cm -2 ) = 21.33 ± 0.12 is quite uncertain. The fit is 
plotted in Fig. 5 (top panel) superimposed to the smoothed 
UVES spectrum. The fit is particularly poor on the wings, 
possibly because more than one component is needed to 
model the absorption. To obtain a better estimate of the 
column density, we used the FORS2 spectrum, that has a 
better S/N (see next section for details). The two component 
fit shown in Fig. 5 (bottom panel) gives a better represen- 
tation of the data. The two components are centered 
at Zi = 1.944 and Z2 = 1.975, respectively, and have 
column densities of log(7V H i/cm" 2 ) = 20.82 ± 0.14 and 
log(iVHi/cm~ 2 ) = 20.79±0.12, respectively. The FORS2 
total column density is log(JV H i/cm" 2 ) = 21.11±0.10. This 
is our best fit result for Nh and will be used in the 
following. The metallicity has been derived summing all 
non saturated component, excited level and ionic contribu- 
tions belonging to the same atom, dividing these values by 
Nh and comparing them to the corresponding solar values 
given in Asplund et al. (2009). The upper limits of compo- 



nent III result in an increment of the total column densities 
of < 20% in the worst cases, so they were not included in 
the computation. The results are listed in Table 4. Column 

2 reports the total abundance of each atom, while columns 

3 and 4 report the absolute and solar-scaled Nx/Nh ratios, 
respectively, with X the corresponding element in column 

I. Lower limits are reported whenever saturation does not 
allow us to securely fit the metallic column densities (see 
e.g., Fig. 4, where the line profiles reach the zero value of 
the normalized flux). In particular, for OI and AM we con- 
sidered also the values of the third, saturated component, 
while for Sill this has not been considered, since the fit 
to the SiIIA1260 resulted in a N value of component III 
which is considerably lower than that of components I and 

II. We derived metallicity values between 0.3 and 0.05 with 
respect to the solar ones. We caution, however, that many 
transitions belonging to other ionization states, which are 
commonly observed in GRB afterglow spectra could not be 
taken into account, because they are outside the UVES- 
dichroic-1 spectral range. In addition, the dust depletion 
can prevent the observation of part of the metallic content 
of the GRB 081008 host galaxy. The reported relative abun- 
dances should then be considered as lower limits to the true 
GRB 081008 metallicity. However, some considerations on 
higher ionization states are possible analyzing the FORS2 
spectra, and dust content can be investigated through the 
study of the depletion pattern (see sects. 5 and 6). 
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Table 3. Absorption line logarithmic column densities for the three components of the main system, derived from the UVES spectrum. 



Species 


Observed transitions 






N (cm" 2 ) 








HI 2 S 1/2 


Lya (UVES) 






21.33 


± 0.12 








HI 2 S 1/2 


Lya (FORS2) 
Components 


20.82 ± 0.14 


(zi = 1.944) 


21.11 ± 0.10 


20.79 ± 0.12 


(» 


= 1.975) 


Metals 


Components 


I (0 km 


b ) 


II (-20 km s -1 ) 


III (-78 km i 




Species 


Observed transitions 


N (cm -2 ) 


b (km s ) 


N (cm -2 ) 


b (km s 1 ) 


N (cm -2 ) 


b (km s -1 ) 


OI 3 P2 


A1302 


> 14.34 


SAT 


> 14.54 


SAT 


> 14.88 




SAT 


A1II 1 S 


A1670 


> 12.82 


SAT 


> 13.26 


SAT 


> 12.78 




SAT 


Aim 2 5i /2 


A1854, A1862 


12.96 ± 0.04 


12 


13.04 ± 0.04 


17 


< 12.3 




N/A 


Sill 2 P° 

1/2 


A1260, A1808 


15.44 ± 0.03 


9 


15.08 ±0.18 


4 


> 13.44 




SAT 


Si IT 2 P° 
Slli ^3/2 






9 


< 15.0 


N/A 


< 15.0 




N/A 


CrII 2 S , 1/ / 2 


A2056, A2062, A2066 


13.51 ± 0.05 


16 


13.55 ±0.04 


21 


< 13.3 




N/A 


tell a"L)g/2 


A2249, A2260 


i rr nn 1 n no 

15.09 ± 0.03 


12 


14.95 ± 0.04 


20 


< 14.7 




N/A 


Fell a 6 £> 7 /2 


A1618, A1621 


13.95 ± 0.05 


12 


< 13.7 


N/A 


< 13.7 




N/A 


Fell a 6 D 5/2 


A1629 


13.71 ± 0.09 


12 


< 13.7 


N/A 


< 13.7 




N/A 


Fell a 6 P> 3/ 2 


A1634, A1636 


13.59 ±0.09 


12 


< 13.5 


N/A 


< 13.5 




N/A 


Fell a 4 P 9/2 


A1612, A1637, A1702 


14.19 ±0.05 


12 


13.76 ±0.07 


20 


< 13.2 




N/A 


Fell a 4 D 7/2 


A1635 


13.35 ± 0.12 


12 


< 13.15 


20 


< 13.15 




N/A 


NiH 2 D 5 /2 


A1741 


13.52 ± 0.08 


12 


13.33 ±0.13 


11 


< 13.3 




N/A 


Nill 4 P 9/2 


A2166, A2217, A2223 


13.57 ±0.02 


12 


13.27 ±0.04 


11 


< 12.9 




N/A 


Znll 2 S 1/2 


A2026, A2062 


12.82 ± 0.06 


4 


12.87 ±0.04 


54 


< 12.4 




N/A 



4.2 Excited levels 

The level structure of an atom or ion is characterized by a 
principal quantum number n, which defines the atomic level, 
and by the spin-orbit coupling (described by the quantum 
number j), which splits these levels into fine structure sub- 
levels. Excited features are routinely detected in GRB ab- 
sorption spectroscopy, at the host redshift, due to the popu- 
lation of both n > 1 and/or n = 1 fine structure levels. GRB 
081008 behaves the same way. In fact, component I features 
the first and second fine structure levels of the Fell ground 
state {a 6 D), the first fine structure level of the Sill 2 P°, and 
the Fell a 4 F 9/2 , Fell a 4 D 7/2 , Nill 4 F 9 / 2 metastable levels 
(the subscript represents the spin-orbit quantum number j). 
Moreover, the Fell a 4 F 7 / 2 and Nill 4 Fg/ 2 excited states are 
also detected in component II (see Table 3 for details). 

There is conspicuous literature on the population of ex- 
cited states in GRB surrounding medium and their detec- 
tion in the afterglow spectra (see e.g. Prochaska, Chen & 
Bloom 2006; Vreeswijk et al. 2007; D'Elia et al. 2010 and 
references therein). There is general consensus that these 
features are produced by indirect UV pumping by the after- 
glow, i.e., through the population of higher levels followed 
by the depopulation into the states responsible for the ab- 
sorption features. This has been demonstrated both by the 



detection of variability of fine structure lines in multi-epoch 
spectroscopy (Vreeswijk et al. 2007; D'Elia et al. 2009a), and 
through the column density ratios of different excited lev- 
els when multiple spectra were not available (Ledoux et al. 
2009; D'Elia et al. 2009b). 

Concerning GRB 081008, the high column density of the 
first metastable level of Fell (a 4 -F 9 /2) with respect to the fine 
structure levels of the ground state can hardly be explained 
with a level population distribution given by a Boltzmann 
function (Vreeswijk et al. 2007), meaning that collisional 
excitations can be safely rejected. The lack of multi-epoch 
spectroscopy does not allow us to completely rule out the 
possibility that the exciting UV flux come from regions of 
high star-formation rates and not from the GRB. In fact, 
fine structure emission lines are present in Lyman-break, 
high redshift galaxies (see Shapley et al. 2003). If we as- 
sume that this flux comes from the GRB, we can estimate 
the GRB/absorber distance, comparing observed column 
densities to those predicted by a time-dependent, photo- 
excitation code for the time when the spectroscopic obser- 
vations were acquired. The photo-excitation code is that 
used by Vreeswijk et al. (2007) and D'Elia et al. (2009a), 
to which we refer the reader for more details. Our equations 
take into account the (47r) _1 correction factor to the flux 
experienced by the absorbing gas described by Vreeswijk 
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Figure 2. The Fell ground and excited absorption features. Solid 
lines represent the two Voigt components, best-fit model. Verti- 
cal lines identify the component velocities. The zero point has 
been arbitrarily placed at the redshift of the redmost component 
(z = 1.9683). g.s. and n* indicate ground state and n-th excited 
transitions, respectively. 



(2011). We assume that the species for which we are run- 
ning the code are at the ground state before the GRB blast 
wave reaches the gas. The GRB flux behavior before the 
UVES observation was estimated using the data in Y10 
(lightcurve and spectral index), with no spectral variation 
assumed during the time interval between the burst and 
our observation. We concentrate on Fell and Sill levels be- 
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Figure 3. The Nill (top panel), A1III and Sill, (middle panel), 
CrII and Znll (bottom panel) absorption features. Solid lines rep- 
resent the two Voigt components, best-fit model. Vertical lines 
identify the component velocities. The zero point has been ar- 
bitrarily placed at the redshift of the redmost component (z = 
1.9683). g.s. and n* indicate ground state and n-th excited tran- 
sitions, respectively. 
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Figure 4. The OIA1302, A1IIA1670, SiIIA1260 ground state and 
SiIIA1264 fine structure transitions. These lines need a three Voigt 
component model to be fitted and are heavily saturated. The zero 
point has been arbitrarily placed at the redshift of the redmost 
component (z = i.9683). g.s. and n* indicate ground state and 
n-th excited transitions, respectively. 




Velocity (km/s) 

Figure 5. The Lya absorption feature at the GRB 08i008 red- 
shift. Top panel shows the single Voigt component, best-fit model 
for the UVES spectrum. Bottom panel shows the double compo- 
nent, best-fit model for the FORS2 spectrum. The UVES fit 
is poor, while the FORS2 one gives a more reliable de- 
scription of Nh- 

cause the Nill ground state has column densities not far 
from the 90% confidence level of \og(N Nil I /cm~ 2 ) — 13.3, 
and thus the uncertainties on such values are high (Table 
3). In addition, Nill ground state is detected only through 
the A1741 transition, because the lower oscillator strengths 
of the A1709 and 1751 lines prevents a detection of these 
features above the 90% level. The initial column densities of 
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Table 4. Metallicity computed from the UVES data. 



Element X 


logiVx/cm 2 


log N x / N H 


[X/H] 


O 


> 15.12 ±0.06 


> -5.99 ±0.13 


> -2.68 ± 0.11 


Al 


> 13.70 ±0.04 


> -7.41 ±0.13 


> -1.86 ± 0.11 


Si 


15.75 ± 0.04 


-5.32 ±0.12 


-0.87 ±0.10 


Cr 


13.83 ± 0.03 


-7.28 ±0.08 


-0.92 ±0.10 


Fe 


15.42 ± 0.04 


-5.69 ±0.13 


-1.19 ±0.11 


Ni 


13.74 ±0.07 


-7.37 ±0.13 


-1.29 ±0.12 


Zn 


13.15 ± 0.04 


-7.96 ±0.13 


-0.52 ±0.11 



the ground states were computed from the observed column 
densities of all the levels of each ion, i.e., we are assuming 
that the species are not excited at t = 0. The initial val- 
ues for Fell and Sill are log(iV Sl /i/cm- 2 ) = 15.63 ± 0.03 
and log(NFeii /cm" 2 ) = 15.21 ± 0.02 for component I, and 
log(N FeII /cm~ 2 ) = 14.98 ±0.04 for component II. Finally, 
the Doppler parameter used as input of this model has been 
left free to vary between 10 and 20 km s _1 , i.e. the range of 
values that best fit the absorption features of components I 
and II. 

Fig. 6 (top) shows the model that best fits the Fell data, 
obtained for a distance of 50 pc and a Doppler parameter of 
20 km s _1 . Fig. 6 (bottom) reproduces the behaviour of the 
reduced \ 2 as a function of the distance GRB/absorber. The 
distance of component I from the GRB explosion site results 
di,Feii = 5ll"iJ pc at the 90% confidence level. The same 
calculation was performed using the Sill atomic data. The 
results are displayed in Fig. 7, and the estimated distance 
is di^sui = 52 ± 6 pc, which is consistent with what was 
estimated using the Fell data. 

For component II, we have much less excited transitions. 
Fig. 8 shows the model which best fits the Fell data and the 
two theoretical curves compatible within the error bars for 
the Fell a 4 Fg/ 2 excited level column density, which is actu- 
ally the only one with a positive detection in component II. 
The resulting distance between the GRB and this absorbing 
component is dn^eii = 200lgQ pc (90% confidence level), 
a larger value than di, as expected given the lack of excited 
transition in component II. 



5 FORS2 SPECTROSCOPY 

In the framework of the ESO program 082.A-0755, we ob- 
served the afterglow of GRB 081008 also with the FORS2 
low resolution spectrograph (R — 780), mounted on 
VLT/UT1. We took three spectra of 900 s each, starting 
around Oct 09 at 00:20 UT (about 4.4 hours after the burst). 
We used the BOOB grism, whose spectral coverage is 330—630 
nm. The extraction of the spectra was performed within the 
MIDAS environment. Wavelength and flux calibration of the 
three spectra were obtained by using the helium-argon lamp 
and observing spectrophotometric stars. Tab. 1 reports a 
summary of our FORS2 observations. We searched for vari- 
ability in the W r of the FORS2 absorption lines, but we 
found none at the 2a level. This is not surprising, since fine 
structure and excited lines are expected to vary by less than 
~ 0.05 decades (in column density) during the acquisition 
time of the FORS2 spectra, which is ~ 15 min rest frame 
(see Figs. 6-8). Since no variability is detected, we co-added 
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Figure 6. Top panel: Fell column densities for the ground level 
(open circle), fine structure levels of the ground state (filled 
circles), first metastable (open square) and second metastable 
(open triangle) transitions for component I in the spectrum 
of GRB 081008. Column density predictions from our time- 
dependent photo-excitation code are also shown. They refer to 
the ground level (dotted line), fine structure level (solid lines), 
first and second excited level (dashed and thick solid lines, re- 
spectively) transitions, in the case of an absorber placed at 50 
pc from the GRB. Bottom panel: the reduced \ 2 as a function of 
the distance for the model reproduced in the upper panel. Dashed 
lines indicate the best fit distance and enclose the 90% confidence 
range. 



15.5 
15.4 

Q2 
I 

| 15-3 

uo 15.2 
o 

15.1 
15 




d=52pc 



80 



90 

trest ( mm ) 



100 



Figure 7. The Sill column densities for the ground level (open 
circle) and first fine structure level (filled circle) transitions for 
component I in the spectrum of GRB 081008. Column density 
predictions from our time-dependent photo-excitation code are 
also shown. They refer to the ground level (dotted line) and first 
fine structure level (thick solid line) transitions, in the case of an 
absorber placed at 52 pc from the GRB. The two thin solid lines 
display the models which enclose the fine structure level data at 
the 90% confidence level (error bars for this transition arc drawn 
both at la and 90% confidence levels). 



the three spectra, to improve the signal-to-noise ratio of our 
data, to obtain a value of ~ 60 - 80 at A > 4000A. The 
resulting spectrum is presented in Fig. 9, together with the 
spectral features identified at z — 1.97, a redshift consistent 
with that estimated using the UVES data. A list of the fea- 
tures detected in the FORS2 spectrum is reported in the 
first column of Table 5. 

The Voigt fitting procedure is not adequate to com- 
pute the column densities of metallic species in low resolu- 
tion spectroscopy. In this case, the Curve of Growth (COG) 
analysis (see e.g. Spitzer 1978) must be applied. For weak 
absorption lines, with width W r < O.lA, and for Doppler pa- 
rameters b > 20 km s _1 , W r is proportional to the column 
density N, and virtually insensitive to the Doppler param- 
eter itself. For stronger lines this does not hold any more, 
and the relation between W r and N is described by a COG, 
which is a function of b. In order to fit the correct COG to 
the data and to estimate b, different transitions (with differ- 
ent oscillator strengths /) of the same species are needed. 
Following Spitzer (1978), we built up a code to perform this 
fit on our FORS2 data. To test our code, we compute W r 
for all the UVES transitions featuring two components, and 
apply our fitting program. The result of the fit is shown in 
Fig. 10 (top panel), and the estimated column densities are 
reported in Table 6 (errors are given at the la level). 

The effective Doppler parameter evaluated from the fit, 
b = 23 km s _1 , is compatible with that estimated using the 
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Figure 9. The flux-calibrated, co-added FORS2 spectrum of the GRB 081008 afterglow, together with the spectral features identified 
at z = 1.9683. 



line fitting profile. To compare the column densities esti- 
mated with the two methods, we sum for each species the 
contribution coming from the two components using the line 
fitting method (see values in Table 3), and report the re- 
sults in Table 6. The agreement between line fitting and 
COG analyses is very good: each column density is within 
la from the corresponding value estimated using the other 
method. The only exception is the FeII5s, whose column 
density values however overlap at the 2a level. 



We now apply the COG analysis to the FORS2 spec- 
trum. First of all, we compute the W r from the data. As 
shown in Table 5, despite the identification of nearly 30 
transitions, reliable W r can be evaluated only for 13 (second 
column of the table). This is because the lower FORS2 reso- 
lution does not enable to separate many of these transitions 
which are blended with each other. We then run the COG 
code using the FORS2 W r , and evaluate the correspond- 
ing column densities. The results are shown in the third 
column of Table 5, while the last column shows the UVES 
column densities for comparison. Errors are again at the la 
level, and Fig. 10 (bottom panel) shows the graphical out- 
put of the fit. The effective Doppler parameter estimated 
(6 = 31 ± 2 km s _1 ) reproduces quite well the combination 
of the values computed for component I and II (~ 10 and 
~ 20 km s _1 , respectively, separated by ~ 20 km s~ x ) using 
the line fitting method. The FORS2 and UVES spectra give 
consistent column densities, with the 3a confidence regions 
overlapping in the worst cases. 



6 CONCLUSIONS AND DISCUSSION 

In this paper we present high and low resolution spec- 
troscopy of the optical afterglow of GRB 081008, observed 
using UVES and FORS2 spectrographs at the VLT ~ 5 hr 
after the trigger. We detect several absorption features (both 
neutral and excited) at the common redshift of z — 1.9683. 
The spectra show that the gas absorbing the GRB afterglow 
light can be described with three components identified in 
this paper as I, II and III, according to their decreasing ve- 
locity values. 

We estimated the distances between the GRB and the 
absorbers. We find a distance for component I of d,Feii,i = 
5lljJ pc and dsm,i — 52 ± 6, using Fell and Sill, respec- 
tively. The Sill leads to a smaller uncertainty because its 
fine structure level is more sensitive to the flux experienced 
by the absorber. Other papers mainly use Fell as distance 
estimator, so for a safer comparison is better to consider 
our Fell value. For component II, this distance is greater, 
du — 200lgQ pc. We stress that these values are obtained 
assuming a three component absorber. However, we can not 
exclude a higher number of components, because our spec- 
trum has a low S/N and a limited resolution. Component II 
is far away from GRB than component I, as expected given 
the lack of fine structure lines in this absorber. Component 
III does not show excited levels at all, and only shows low 
ionization states. Therefore, this is produced by an absorber 
located even farther from the GRB, in a region which is not 
significantly influenced by the prompt/afterglow emission. 

Component I of GRB 081008 is the closest to a GRB 
ever recorded. In fact, for the 6 other GRBs for which the 
GRB/absorber distance have been estimated, the closest 
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Figure 8. The Fell column densities for the ground level (open 
circle), first fine structure level (upper limit) and first excited 
level (open square) transitions for component II in the spec- 
trum of GRB 081008. Column density predictions from our time- 
dependent photo-excitation code are also shown. They refer to 
the ground level (dotted line), first fine structure level (solid line) 
and first excited level (thick dashed line) transitions, in the case 
of an absorber placed at 200 pc from the GRB. The two thin 
dashed lines display the models which enclose the excited level 
data at the 90% confidence level (error bars for this transition 
are drawn both at la and 90% confidence levels). 



components are at d = 80 — 700 pc from the GRB (Vreeswijk 
et al. 2007; D'Elia et al. 2009, a,b; Ledoux et al. 2009; D'Elia 
et al. 2010; Thone et al. 2011). The values reported in liter- 
ature have been corrected for the 4(7r) _1 ^ 2 factor discussed 
by Vreeswijk (2011). This behaviour can be interpreted as 
due to a dense environment close to the GRB explosion site. 
This high density is possibly witnessed by the a non neg- 
ligible dust amount (see below) and by the metal content 
of the GRB surrounding medium. In fact, the GRB 081008 
surroundings have the highest metallicity and the highest 
abundances of, e.g., Fell and Nill, among this sub-sample 
of GRBs. This high density in the GRB surroundings could 
constitute a barrier to the GRB prompt /afterglow emission, 
that is not able to strongly excite the interstellar medium 
up to the distances reached by the other GRBs. 

The neutral hydrogen column density is 
log(A?H,opt/cm~ 2 ) = 21.11 ± 0.10, while that estimated 
from Swift XRT data is log Af H ,x/cm" 2 = 21.66to.26 
(Campana et al. 2010). The latter value is for a solar abun- 
dance medium. Using A^H, op t we evaluate the GRB 081008 
host galaxy's metallicity. The values we find are in the 
range [X/H] = —1.29 to —0.52 with respect to the solar 
abundances. This value lies in the middle of the GRB dis- 
tribution, (Savaglio 2006; Prochaska et al. 2007; Savaglio, 
Glazebrook & Le Borgne, 2009). From X-ray data a limit of 
[X/H] > —1.83 (90% confidence limit) can be set assuming 



Table 5. GRB081008 absorption features detected in the FORS2 
spectrum, together with their W r and column densities. UVES 
data are shown for comparison. 



Transition 


w r (A) a 


N b 


UVES N 6 


CIIA1334 


BLEND 


- 


- 


CIVA1548 


BLEND 


- 


- 


CIVA1550 


BLEND 




- 


OIA1302 


BLEND 




SAT 


OIA1304 


T"»T T71TVTT" \ 

BLEND 






OIA1306 


T"»T nMFl 

BLEND 






A1IIA1670 


0.74 


1 4 60 


SAT 


A1IIIA1854 


0.30 


13.42 ± 0.02 


13.30 ± 0.03 


A1IIIA1862 


0.17 


13.42 ± 0.02 


13.30 ±0.03 


SiIIA1260 


BLEND 


15.74t°, ; ii 


15.60 ±0.04 


SiIIA1304 


BLEND 


15.74±°;i; 


15.60 ± 0.04 


SiIIA1526 


0.67 


15 74+°- 15 


15.60 ±0.04 


SiIIA1808 


0.25 


15 74+0.15 


15.60 ± 0.04 


Sill A1264 


BLEND 


15.17 ±0.01 


15.21 ± 0.05 


Sill A1309 


BLEND 


15.17 ±0.01 


15.21 ± 0.05 


Sill A1533 


BLEND 


15.17 ±0.01 


15.21 ± 0.05 


Sill A1816 


0.07 


15.17 ±0.01 


15.21 ± 0.05 


SilV A1393 


0.47 


- 


- 


SilV A1402 


0.40 


- 


- 


CrIIA2056 


0.27 


13.94 ±0.02 


13.83 ± 0.03 


CrIIA2062 


BLEND 


13.94 ±0.02 


13.83 ± 0.03 


CrIIA2066 


0.17 


13.94 ±0.02 


13.83 ±0.03 


FeIIA1608 


0.54 


15 29+0' 13 


15.33 ±0.02 


FeIIA1621 


BLEND 




14.95 ± 0.04 


FeIIA1629 


BLEND 




14.95 ± 0.04 


FeIIA1702 


0.24 


14.16 ±0.02 


14.33 ±0.05 


ZnIIA2026 


0.24 


13.22 ±0.02 


13.15 ± 0.04 


ZnIIA2062 


BLEND 


13.22 ± 0.02 


13.15 ± 0.04 



a The error on W r is O.OlA(ler). 

6 All values of the column densities are logarithmic (in cm -2 ). 

a solar abundance pattern and requiring that the absorbing 
medium is not Thomson thick. If we set the metallicity to 
[X/H] = —0.5, the absorbing column density in the X-rays 
is higher, namely, log Nu,x/cm~ 2 = 22.24 + . ' 30 and higher 
for lower metallicities. Fynbo et al. (2009) and Campana 
et al. (2010) show that in GRBs with a detectable Lya 
feature (i.e., those at z > 2) Nh,x is on average a factor 
of 10 higher than N H , op t, and GRB 081008 follows this 
trend. The intense GRB flux, which ionizes the hydrogen 
and prevents part of it to be optically detected, is the 
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Table 6. Comparison between UVES column densities cvalutated 
with the line fitting and COG methods. 



7.6 7.8 8 8.2 

Log NfX (cm- 1 ) 



-3.5 



\ 

u 



o 



-4.5 -- 



n — i — i — | — i — i — i — i — | — i — rp — i — | — i — i — i — r 

b = 31 km/s / 



8.4 




8.5 9 9.5 

Log NfX (cm- 1 ) 



Figure 10. Top panel: the COG analysis tested using the UVES 
species featuring two components. Bottom panel: COG analysis 
applied to the FORS2 lines with a measured W r - Solid lines repre- 
sent the best fit obtained using the reported b values. Dashed lines 
show the b — oo curve for comparison. The COG fits component 
I and II together in both plots. 



common explanation for this discrepancy (Fynbo et al. 
2009; Campana et al. 2010; Schady et al. 2011). 

It is worth noting that observed abundances of Fell and 
Znll are significantly different ([Fe/H]= -1.29 ± 0.11 and 
[Zn/H]= -0.52 ± 0.11). This can be ascribed to the differ- 
ent refractory properties of the two elements, with the for- 
mer that preferentially tends to produce dust grains while 



Specie 


N (COG analysis) 


N (Line fitting) 


A1III 


13.29tg-°g 


13.30 ±0.03 


Sill 


15.661°;°° 


15.60 ±0.03 


CrII 


1 9 QQ+0.03 
±o.oo_ 07 


13.83 ±0.03 


Fell (g.s.) 


15 S1+°' ()l 

io.oi_ 07 


15.33 ±0.02 


Fell a 4 F g/2 




14.33 ±0.05 


Nill (g.s.) 


l3.8l±g;gS 


13.74 ±0.07 


Nill a 4 F 9/2 


,„ 7 „+0.05 
io.(o_ og 


13.75 ±0.02 


Znll 


,„ , ,+0.05 


13.15 ±0.04 



All values are logarithmic (in cm 2 ). 



the latter prefers the gas phase. The comparison between 
these 'opposite' elements can thus provide information on 
the dust content in the GRB environments. In order to 
be more quantitative, we derive the dust depletion pattern 
for the GRB 081008 environment, following the method de- 
scribed in Savaglio (2000). We consider the four depletion 
patterns observed in the Milky Way, namely, those in the 
warm halo (WH), warm disk + halo (WHD), warm disk 
(WD) and cool disk (CD) clouds (Savage & Sembach 1996). 
We find that the best fit to our data is given by the WH 
cloud pattern, with a metallicity of log Zgrb /Zq ~ —0.5 
and a GRB dust-to-metal ratio comparable to that of the 
WH environment, e.g., d/dwH = 1 (Fig. 11). This metallic- 
ity value is consistent with our [Zn/H] measurement. This 
agreement is self-consistent with the use of zinc as a good 
indicator of metallicity. Since the latter quantity is linked 
to the extinction (see e.g., Savaglio, Fall & Fiore 2003) we 
derive A v ~ 0.19 mag along the GRB 081008 line of sight. 
We check this value by modeling the flux-calibrated FORS2 
spectrum. The SED is dominated by the Lye* which is dif- 
ficult to model given the high fluctuations and other ab- 
sorption lines. Anyway, the inferred Av value is low and 
compatible with that evaluated from the dust depletion. 

Another hint of dust is the non detection of Fell in the 
third component. The Fell column densities in components 
I and II are very similar, and this lets us believe that in 
component III the iron is present as well, but in the dust 
form. The higher presence of dust in components far away 
from the GRB has already been pointed out by D'Elia et al. 
(2007) . They report a possible presence of dust in component 
III of GRB 050730, while the closer component II (featuring 
Fell fine structure lines) shows more iron in the gas state. 
The detection of more dust far away from the GRB can 
be explained since dust grains containing iron tends to be 
efficiently destroyed during a blast wave occurring after a 
GRB explosion (Perna, Lazzati & Fiore 2003). 

The analysis of the FORS2 spectra extends our sur- 
veyed wavelength range, allowing the detection of higher ion- 
ization species, such as CIV and SilV. Anyway, line profiles 
of high and low ionization species rarely match in redshift 
space and often if they do, it is because the line blending 
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Figure 11. Depletion patterns in the absorbing gas of 
GRB 081008. Filled squares are taken from average gas-phase 
abundance measurements in warm halo (blue), warm disk + halo 
(green), warm disk (red) and cool disk (cyan) clouds of the Milky 
Way (Savage & Sembach 1996). Filled circles represent our data 
points, which are best fitted by the warm halo cloud pattern. 

cannot be resolved in a spectrum, regardless of resolution 
and S/N. 

We stress that the availability of simultaneous high and 
low resolution spectra of a GRB afterglow is an extremely 
rare event. In this context, the comparison of the column 
densities obtained fitting the line profile of a high resolu- 
tion spectrum with that estimated by the Curve of Growth 
analysis applied to a low resolution one could be extremely 
important. In fact, this can help to determine a range of 
column densities for which it is safe to apply the Curve of 
Growth analysis when high resolution data are missing. This 
is because high column densities can result in the saturation 
effect, a problem that is difficult to address using low resolu- 
tion spectra only (see e.g. Penprase et al. 2010). Prochaska 
(2006) widely discuss the limits and perils of the COG analy- 
sis applied to low resolution data. They find that this kind of 
analysis tends to underestimate the column densities of the 
absorbing species. This is because strong transitions drive 
the COG fit since the relative error associated to their W r is 
smaller than that for weak ones. Nevertheless, strong transi- 
tions are more affected by saturation, and in order to match 
their observed column densities, the COG fit is forced to- 
wards high values of the effective Doppler parameter. High 
resolution data often show that the main contribution to the 
column density of strong transitions comes from one narrow 
component. On the other hand, the main contribution to the 
W r comes from other components which account for a small 
fraction of the column density. These inferred high values for 
the effective Doppler parameter are thus mimicking a more 
complex situation, with the result of underestimating the 
real column densities. For what concerns GRB 081008 the 
UVES observations show no or just mild saturation even for 
the strongest transitions, and the two main components give 
a similar contribution to the total column densities. This is 
the reason why there is a good agreement between COG 



analysis of low resolution data and line fitting analysis of 
high resolution ones for this particular GRB (within 3<r in 
the worst cases). 

Finally, we detect two weak intervening systems in our 
spectra. The first one is a CIV absorber in the FORS2 spec- 
trum at z = 1.78, and the second one is a Mgll system 
in the UVES spectrum at z = 1.286. This last system has 
W r (MgII\2796) = 0.3A, the detection limit being O.lA at 
the 2cr confidence level. The redshift path analyzed for Mgll 
is z = 0.18 - 0.38 and z = 0.71 - 1.43 for the UVES spec- 
trum, and z = 0.36 - 1.21 for the FORS2 one. 
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